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Abstract 

r-i 

In reference [l[, Michel Gaudin showed the Wick's theorem at finite temperature using a di- 
agonal Hamiltonian. We extend the Gaudin's prove for a statistical density operator which de- 
pend on a quadratic Hamiltonian. To illustrate the utility of the theorem, we evaluate the ratio 
J [{N 2 ) — (N) 2 ]/(N) of a homogeneous weakly interacting Bose gas at temperature below the Bose- 

Ft" 1 

Einstein condensation temperature. At this condition, the quadratic Hamiltonian approximation 
is valued and, in this evaluation, we show the sub-Poissonian behaviour of the fundamental state 
distribution at zero temperature. 
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I. INTRODUCTION 



We are interested in describing a quantum system enclosed in a volume V liable to 
the energy and particle reservoir. Such system branches of science as quantum gases and 
liquid, nuclear models, quantum information and cosmological dark matter. Attention to 
these systems is primarily due to the possibility of a deeper understanding of quantum 
thermodynamics in open systems. 

At zero temperature, the system can be represented by a symmetric (regarding to boson 
particle) or antissymmetric (regarding to fermion particles) state vector, \^(t)), in the Fock 
space. The dynamics of the system is achieved by the solution of the Schrodinger equation, 
\^(t)} = e~ ll £~\fy}, where H is the Hamiltonian operator and \^} is the initial state vector. 

In contact to the energy and particle reservoir, the system has a finite temperature and 
can be described by the statistical density operator in the Fock space, T>, whose properties 
are hermitian, positive defined and the trace equal to unity. Hence, the thermodynamics 
variables are determined by the ensemble average (O) where O is the quantum observable 
related to the thermodynamics variable and is obtained by the prescription, 

(6) = TrT>6 (1) 

where the trace is over the Fock space. 

In statistical mechanics, the statistical density operator is determined by the principle 
of maximum entropy which is defined in terms of the statistical operator as the ensemble 

average, 

S=-k B TrV\nV (2) 

where the k B is the Boltzmann constant. Whereas the system is in contact to the energy 
and the particle reservoir, we have introduced Lagrange multipliers, T corresponding to 
the temperature and ji corresponding to the chemical potential, to constraint the ensemble 
average energy E and the particle number average N. Hence, we are led to minimize the 
functional, 

4>{V) = TrV(k B T \nV + )C) (3) 

where is the functional called the thermodynamic potential and K, — H — jj,N is the 
canonical Hamiltonian with N as the particle number operator. The variational equation 
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for V yields the statistical operator as, 



e -pK 



V = =r. (4) 

Tre-f> K 

where (3 = l/ksT. 

Without loss of generality and for the sake of the simplicity to take the trace in the 
ensemble average, we build any operators, O, as linear combination of the products of the 
creation and annihilation operators a\ and «j of the state vectors in the Fock space. Thus 
we can written O as, 

Q = ^2 Unm{®Wj ■ ■ • a\pii OiiOLj ■ • ■ a k ai ), (5) 

n,m s v ' * 

n m 

where U nm is a constant. The operators a\ and a>i satisfy the commutation rules, 

[cKj, aj] e = CKjO;] — ea^oti = (6) 

where e equal to +1 when a\ creates a boson state and —1 when a, creates a fermion state. 

In order to proceed the trace in the average prescription, it is appropriated to approach 
the canonical Hamiltonian to the quadratic Hamiltonian 

fc*ff = 2 X^M n ' n / + a i a \) + A *j a i a j + Arafat], (7) 

where the elements hij and Ay are dynamics parameters. 

In the quadratic form, the Hamiltonian can be diagonalised by a canonical Bogoliubov 
transformation in creation and annihilation operators, a\ and ctj, 

t — v t _ 

where the sign j consist of time-reversed states and the transformation parameters, Uij and 
Vij, obey the normalization conditions, 

E(h*l 2 -eM 2 ) = i, ( 9 ) 

ensuring the transformation be canonical. Such parameters are given to diagonalize the 
canonical Hamiltonian and be possible to write the quadratic Hamiltonian, (jHJ), in a diagonal 
form, 

t eff = J2 E ivh- (10) 



The generalized Wick's theorem demostrated in [l| asserts the ensemble averages is equal 
to the sum over all possible contracted terms, 

(O) = ^2U nm ((a\a] ■ • ■ ot\a\oaotj ■ ■ ■ a k ai}) 

n,m 

= U nm {a\a)(a\a ] m ■ ■ ■ a k a t ) + ea\a\(a)a ] m ■ ■ ■ a k ai) H ) (11) 

n,m 

where the contracted terms are given by, 



ajaj = 



(12) 



Then the ensemble averages can be written as linear combination of ensemble averages of 
the product of one creation and other annihilation operators. 

In the section [Til we demonstrate the generalized Wick's theorem at finite temperature 
using the quadratic Hamiltonian (JTJ). Differently from the generalized Wick's theorem de- 
mostrated in |l|, the contracted terms depend on the transformation parameters Uij and Vij. 
In section II III we use this theorem to show the sub-Poissonian distribution in the particle 
number ensemble average to a homogeneous weakly interacting Bose gas. We conclude the 
subject in section [TV] 

II. DEMONSTRATION OF THE WICK'S THEOREM 

Let a ensemble average be {abc ■ ■ ■ ef) where a, b and etc. are creation or anihilation 
operators. The generalized Wick's theorem asserts, 

(abc ■ ■ ■ ef) = ab(bc ■ ■ ■ ef) + eac(bd ■ ■ ■ ef) + • • • . (13) 

To proceed the demonstration, we assume the operator a as a\ and write the trace in 
terms of transformed operators, 77} and m, 



Tr{ajbc ■ ■ ■ eft)} = £ ^Tr^jfec • ■ ■ efV} - ^Tr{^foc ■ • • efV}. (14) 
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whereas, 77]© = e l3Ej 'Dri^ and rjj'D = e l3E: >'Dr]j, due to the canonical hamiltonian, lC e ff, 
being diagonal in this news operators, the trace depend on, 



Tr{6c • • • efrj}V} = e^ E 'Tr{r]}bcd ■ ■ ■ efV}, (15) 



Then we have 



Ti{q}bcd ■ ■ ■ efV} 

jS^k) Tl{cd ■ ■ ■ efi)} + e jr^k) Tr{M • • • efi)} + 



and 



Trfabcd ■ ■ ■ efV} = 

[% '' b ^ -Ti{cd ■ ■ ■ efV} + e - — °} 6 Ti{bd ■ ■ ■ eft)} + 



Using (jT6]) and (jTTj) in (jl4j) . we have 
Ti{a\bc ■ ■ ■ efV} = 



= E ( ^MA ; _ vulval ) Tr{cd ... efi)} + 

+ E f 7T^#V ~ ^ {\ ^ {bcd ■■■eV}, 



in the same way using a = oij, we have 



Tr{aj6c • • • e/X>} 

Vf W fr&jK |Tr{cd---efP| + 

+ e E(, [ ^ C] ^L-^ t,C] ^NM---e/P} + - 
^ l (l-ee-^) (1-eeW' 1 1 



^ 1 (1 - ee"^ ) (1 - ee^ ) / 1 1 



Then, we show the contracted terms is 



and, therefore, we write the ensemble average as 

Tr{abc ■ ■ ■ efV} = abTr{cd ■ ■ ■ efV} + edcTr{bd ■ ■ ■ efV} + 

+ (^Tr{6c • • • efV} + ■ ■ ■ + afTr{bcd ■ ■ ■ eV} (21) 

and, thus, we demonstrated the generalized Wick's theorem at finite temperature using the 
quadratic Hamiltonian whose the contracted terms are given in ( 120]) and depend on the 
transformation parameters. 

III. AN ILLUSTRATIVE THEOREM APPLICATION 

To illustrate the utility of the generalized Wick's theorem, we evaluate, in the homoge- 
neous weakly interacting Bose gas at low temperature, the parameter Q define by, 



(22) 

which is the measure to the departure of the particle number distribution from a Poissonian 



distribution (4]. 

For the ideal Bose gas, the fluctuation of the particle number, at temperature below the 
the Bose-Einstein condensate transition, have been studied considering the gas inside a box 
[5I or in harmonic trap jfj. In both, it is shown that microcanonical fluctuations below 



the critical temperature tend to zero and scale with the number of particles as l/y(N). 
However, in the case of the weakly interacting gas, it is necessary consider the Hamiltonian 
with interacting terms. 

We evaluate the ensemble average of product of four boson operators, a, b, c e d. Using 
(f2~TT) we have 



Tr {abcdV} = abTr{cdV} + dcTr{bdV} + adTr{bcV}, (23) 

whereas 

Tr {abcdV} = Ti{abV}Ti{cdV} + Tr{acV}Tr{bdV} + Tr{adV}Tr{bcV}, (24) 

we show the expansion of the four operator ensemble average in the sum of the two operator 
ensemble average through the Wick's theorem. 
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Thus the evaluation of the parameter Q at low temperature regime depends on the fluc- 
tuation ensemble average, (N 2 ), and the particle number ensemble average, (N). Expanding 
the creation and annihilation operator in a basis of the plane wave functions with momentum 
k, the fluctuation is, 

(N 2 ) =J2Tr{ala n al,a n ,V}, (25) 

k,k' 

and the particle number ensemble average is, 

(A0 = £Tr{4a fe ^}. (26 ) 
k 

When the temperature is below to transition, the Hamiltonian can be approached to a 
quadratic and the vacuum can be represented by a coherent state characterized by complex 
number zq. This vacuum can be described by a shift in particle operator 

a k = c k + z ok,o 



(I 



k °k 



4 + z ok, 



(27) 



then, using the generalized Wick's theorem, we determine the average values of shifted 
operators. 

The fluctuation can be written as, 



<^ 2 > = E<4 C ^> +2 ^E<* S > + 

k,k' k 

+ ^o((44) + (co4) + (eta) + (coco)) + 4 (28) 
Using the generalized Wick's theorem, we obtain 

k,k' 

+ 2z 2 Q E(4 c fc) + ^o(l + 2(cta> + (44) + (coco)) + 4 (29) 

Taking the Bogoliubov transformation, 

eg = cosh((7£)&£ - sinh(ag)6 f _. (30) 

where the hyper-geometric angle is adjusted to diagonalize the quadratic Hamiltonian 
(J7J). We have the contracted terms as function of the hyper-geometric angle, 



i C k C k') ~ 2 



coth (£ft) (1 + sinh(2ag)) 
coth (^f-) sinh(2cr^) - 



(31) 



where Ef. is the low-lying excitation energy. 

Hence, the fluctuation can be written as function of tthe excitation energy, temperature 
and hyper-geometric angle, 



(TV 2 ) 



+ ~£[coth( 

k \ 



0E k 



sinh(2(7?) - 1] + 



k 

+ Zq lim 

fc^O 



coth 
coth 
coth 



(l + sinh(2^)) 
sinh(2cxg) — 1 
sinh(2cr^) 



+ 



+ 



(32) 



In the cases of the fluctuation convergence and at zero temperature the expression ( 132|) 
resume to, 



1 



(33) 



Considering cr = and knowing that (sinho^) 2 < ( smh2gfc ) we can define a minimum 
limit to (N 2 ) and, therefore 

1 



Q > 



(N) 



(34) 



Although the vacuum is a coherent state, the fundamental state show to be narrower 
than the Poissonian distribution due to the Q , for high particle number, goes to zero slower 
than the Poisson distribution. 



IV. CONCLUSION 

We demonstrated the generalized Wick's theorem at finite temperature considering a 
quadratic Hamiltonian and we illustrated this result with a application to evaluate the 
fluctuation particle number of a weakly interacting Bose gas below the condensate transition 
temperature. 

Different from the reference 1|, we prove the generalized Wick's theorem for a quadratic 
Hamiltonian is similar to the simplified Gaudin's demonstration. However, the contracted 
terms expressions is a linear combination of to the Bogoliubov parameters weighted by the 
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commutation rule and the statistical distribution which depends on the symmetry of the 
Fock states. 

An illustrative example was done for determining the parameter Q measure how narrow is 
the particle distribution. The result showed the sub-Poissonian character of the distribution 
due to the Bogoliubov parameters. An appropriate explanation for this behavior will be 
shown in a future work. 
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